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We show that models for homogeneous and heterogeneous nucleation of D-dimensional droplets 
in a d-dimensional medium are described in mean-field by a modified Smoluchowski equation for 
the distribution N(s, t) of droplets masses s, with additional terms accounting for exogenous growth 
from vapor absorption, and injection of small droplets when the model allows renucleation. The 
corresponding collision kernel is derived in both cases. For a generic collision kernel K, the equation 
describes a clustering process with clusters of mass s growing between collision with s oc and 
injection of monomers at a rate I. General properties of this equation are studied. The gel criterion is 
determined. Without injection, exact solutions are found with a constant kernel, exhibiting unusual 
scaling behavior. For a general kernel, under the scaling assumption N(s,t) ~ Y(i) _1 f(s/S(t)), we 
determine the asymptotics of S(t) and Y(t), and derive the scaling equation. Depending on (3 and 
K, a great diversity of behaviors is found. For constant injection, there is an asymptotic steady 
state with N(s, t = oo) oc s~ T and r is determined. The case of a constant mass injection rate is 
related to homogeneous nucleation and is studied. Finally, we show how these results shed some 
new light on heterogeneous nucleation with d = D. For d = D — 2 (discs on a plane), numerical 
simulations are performed, in good agreement with the mean-field results. 



I. INTRODUCTION 

Some very simple and practically important physical 
phenomena can be related to aggregation models, and 
as a consequence a host of experimental, numerical and 
theoretical studies are to be found in the literature. Clas- 
sical fields of application are atmosphere sciences, mate- 
rial sciences, chemical engineering and cosmology, among 
others fflj^]- One of these phenomena is dropwise con- 
densation on a substrate || , for instance water on a cold 
window pain, which bears on important implications in 
heat transfer engineering and material sciences, and gen- 
erates fascinating droplets patterns, also called breath 
figures H. Droplets grow from vapor, and when two 
droplets come into contact, they coalesce to form a big- 
ger droplet, with mass (or volume) conservation. The 
underlying physics is rich, and it is not our purpose to 
fully discuss it here (see ||). Simple models have been 
introduced to describe the late stage of droplet growth 
and coalescence, of basically two kinds |]^,f7|-|l"T| . On the 
one hand, in heterogeneous nucleation models, one starts 
from a fixed number of nucleation sites (in physical situ- 
ations these might be dust particles, substrate defects...). 
Droplets grow on these sites through vapor absorption, 
and when two neighboring droplets overlap, they coa- 
lesce to form a single droplet, thus reducing the number 
of droplets. On the other hand, in homogeneous growth 
models fjj], nucleation can occur anywhere on the sub- 
strate: some very small droplets are randomly deposited, 
which leads to the growth of existing droplets, and the 
creation of new small droplets if deposition occurs in a 
free zone. 

Both kinds of models are aggregation models, with 



the unusual feature that the aggregating particles grow 
between coalescence events. Homogeneous growth has 
the additional feature that there is a source of small 
droplets, which relates this model to aggregation in pres- 
ence of a source, which is also widely-studied in the lit- 
erature |l^-|2l]]. Experiments and numerical simulations 
HU snow that the time dependent droplet mass distri- 
bution N(s,t) (s being the volume of the droplets), ex- 
hibits dynamic scaling, i.e. that for large times N(s, t) oc 
S(t)~ e f(s/S(t)). S(t) is a typical droplet mass (pro- 
portional to<s 2 >/<s> 2 ) and has a power law 
divergence at large time S(t) cx t z . 6 and z are dy- 
namic exponents that do not depend on the fine de- 
tails of the model but only on its main features, such 
as its conservation laws. In heterogeneous growth mod- 
els, the scaling function f(x) is narrow-shaped, whereas 
in homogeneous growth one observes the superposition 
of a monodispersed distribution of large droplets (with 
masses of the same order as S(t)), and a polydispersed 
distribution of small droplets with a power law diver- 
gence of the scaling function f(x) oc x~ T at small x. r is 
nontrivial and less than 6 Q . Such nontrivial polydisper- 
sity exponents frequently appear in aggregation models. 
For usual models (without injection or growth), they can 
be accounted for by Smoluchowski's mean-field approach 
p2j: neglecting fluctuations and multiple collisions, one 
can write down a rate equation, 

d t N(s, t) = ±f N(si,t)N(s - s 1 ,t)K(s 1 ,s - si) d Sl 

- N(s, t) J AT(si,t)if(s, si) ds u (1.1) 

where the collision kernel K(si,S2) is the probability of 
a coalescence event between two droplets of masses S\ 
and S2- Smoluchowski's approach is valid above an up- 
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per critical dimension which is often 2, but is in princi- 
ple model-dependent p3]| . Van Dongen and Ernst 
classified the kernels according to their homogeneity and 
asymptotic behavior: 



suggest that including pair correlations may be necessary 
to find a consistent scaling for homogeneous growth. 



K(bx,by) = b x K{x,y), 
K{x,y)^x^y u {y>x). 



(1.2) 
(1.3) 



Nontrivial polydispersity exponents appear in the case 
[i = J24] ^(| , whereas for > 0, r is equal to 1 + A, and 
for /x < the distribution is bell-shaped. 

In this article we are concerned with studying the 
mean-field Smoluchowski approach to those aggregation 
processes, such as droplet nucleation, where individ- 
ual particles or clusters grow between collisions (with a 
growth law s cx s' 3 ), and/or monomers are injected with 
a possibly time-dependent injection rate. Our first moti- 
vation to this study is a better understanding of droplets 
deposition, growth and coalescence models, and through- 
out the article, the example of both homogeneous and 
heterogeneous droplet nucleation will be used as an il- 
lustration of the results. Certain results presented here 
have already appeared in a summarized form in ]27| ]. 

Sec. [il] discusses nucleation models introduced by Fam- 
ily and Meakin Q , and derives the corresponding Smolu- 
chowski equation under the mean-field assumption. This 
leads us to a generalized Smoluchowski equation with an 
additional exogenous growth term d s (s /3 N) in the left 



hand side of Eq. (1.1) and a time dependent source term 
I(t)S(s — so) in its right hand side. 

Sec. Ill is a general study of the extended Smolu- 



chowski equation corresponding to aggregation with ex- 
ogenousgrowth and injection, such as the one obtained 
in Sec. O, with a generic homogeneous kernel. The gela- 
tion crite rion is investigated, depending on A defined in 
Eq. ( |l.2|) . It is found that the system is nongclling for 
max(A,/3) < 1. Some exact solutions are found in the 
absence of injection (K = 1, with f3 = or j3 = 1). 
The scaling properties of the equation without injection 
are investigated. It is found that the scaling behavior 
depends on A and fi, and is richer than for standard 
Smoluchowski's equation. The occurrence of polydisper- 
sity exponents is discussed. When polydispersity occurs, 
the scaling equation is the same as for standard Smolu- 
chowski's equation, and the methods recently introduced 
by the present authors p6| can be used to compute non- 
trivial r exponents. 

For constant injection of monomer, it is found that 
the distribution reaches at infinite time a polydispersed 
steady state with a power law large s decay, with r = 
(3 + A)/2 if 13 < (1 + A)/2, and r = 2 + A - (3, if 
f3 > (1 + A)/2. Then we consider a constant mass injec- 
tion rate, with a self-consistent I(t), for A = 2(3—1, which 
is relevant to homogeneous nucleation, and we show that 
the injection rate I{t) is vanishing, in agreement with the 
droplets model. We also investigate scaling solutions, and 



Sec. IV is an application of the scaling mean-field re- 
sults to heterogeneous growth with d = D. Droplets 
radii r = s 1 / grow as r « r" (P = 1 + — 1)/D), 
and polydispersity with a nontrivial r occurs for oj > 0, 
while the scaling function is monodipersed for u) < 0. 
Mean-field polydispersity exponents are computed using 
the variational method introduced in |2q|, and described 



in Sec. [II. Numerical results for the scaling function arc 
in qualitative agreement with mean-field results, and the 
expected cross-over from monodispersity to polydisper- 
sity at lo = is observed. 



II. DROPLETS DEPOSITION, GROWTH AND 
COALESCENCE IN MEAN-FIELD 



As mentioned in the introduction, interest in droplet 
nucleation computer models was primarily aroused by 
practical applications in heat transfer engineering (see 
references in [^)). In the last ten year, however, and 
since the seminal work of Beysens and Knobler || , the 
focus was set on the formation of breath figures (see Fig. 
|j|), with computer models aimed to study the kinetics 
of the droplet mass distribution the asymptotic 

surface (or line) coverage [ [Tof , or the time evolution of 
the "dry" fraction (the surface fraction which has never 
been covered by any droplet) JyJ. 

In this article, we shall consider the specific models 
introduced by Family and Meakin |7|||], both for homo- 
geneous and heterogeneous nucleation. We shall now de- 
scribe these models and derive the corresponding Smolu- 
chowski equation. These equations are special cases of a 
general ized Smoluchowski equation which will be studied 
in Sec. lill. 



A. Homogeneous nucleation 

First, let us describe homogeneous nucleation, or the 
deposition and coalescence model. Between t and t + St 
a small droplet of mass sq is randomly deposited on the 
c?-dimensional substrate where it forms a spherical cap 
with radius Sq D '. If it overlaps an existing droplet of 
mass s, they coalesce to form a new droplet with mass 
s + s and radius (s + so) 1 ^ , centered at the center of 
mass of the two coalescing droplets. If the new droplet 
overlaps a surrounding droplet, they coalesce with the 
same rule, and so on. The distribution of droplet masses 
N(s,t) exhibits dynamic scaling, 



N(s,t)~S(t)- e f(s/S(t)). 
The typical mass scale S(t) can be defined by, 



(2.1) 



2 



S(t) = 



< s > 



Js 2 N(s,t) 
fsN(s,t) ■ 



(2.2) 



The dynamical exponents 9 and z are easily determined 
from physical arguments ]5|,|| . 

Since the mass injection rate is constant and the mass 
is conserved in the coalescence process, we must have, 



t oc / sN(s,t)dt oc S(ty 
Jo 



+00 



xf(x)dx, (2.3) 
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which, from the definition of z, implies the scaling law 
z (2 — 9) = 1. Then we note that the fraction of substrate 
"area" occupied by the droplets is, 



FIG. 1. S(t) from numerical simulations of deposition and 
coalescence in d = 1. During one time step, 10 7 droplets of 
radius 0.75 were deposited on a line of length 5 x 10 5 . Ten 
samples were averaged to obtain these data. The scaling laws 
S(t) z with z — D/(D — 1) are very well obeyed. 



r + CG 

/ s*N{s,t)dsoz S{t) 1+ i~ e , (2.4) 
Jo 



and cannot diverge or vanish, so that 9 = 1 + d/D. From 
the scaling law, we get z — D/(D — d). 

The scaling behavior of the total number of droplets 
n(t) depends on the small x behavior of the scaling func- 
tion f(x). If the scaling function is integrable in zero, it 
is easily seen that n(t) oc S 1 ^) 1-9 , whereas if f(x) oc x~ T 
with r > I, n(t) oc S(t) T - e . 

Numerical results as obtained in Q confirm the scal- 
ing hypothesis, and the theoretical values of 9 and z. 
Actually, authors found some values of z a few percents 
smaller than D/(D — d), since simulations did not reach 
sufficiently large times. To illustrate our discussion, we 
present here our own numerical simulations in one dimen- 
sion, for D — 2,3,4. Droplets of radius r$ = 0.75 were 
randomly deposited on a one dimensional "substrate" of 
length L = 5 x I0 5 . During one simulation time step, 
about 10 droplets were deposited. Ten samples were 
averaged to obtain the final data. Fig. |l| shows the evo- 
lution of S(t). We observe a power law S(t) oc t z with 
numerical values of z equal to their theoretical values 
with excellent accuracy (for instance for D = 3 we find 
z = 1.4995±0. 00015 to be compared with the theoretical 
z = 3/2). 
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FIG. 2. Scaling of the mass distribution N(s,t) for droplet 
deposition with d — 1 and D = 3. The picture shows the 
excellent data collapse, with the theoretical value 4/3 for 8, 
of the distribution at four different times when S has re- 
spectively reached the values 6053, 11116, 15539, and 17112. 
The scaling function is composed of a polydispersed contri- 
bution of small droplets and a monodispersed contribution of 
droplets of mass of order S(t). 

The scaling for the distribution function N(s, t) is illus- 
trated in Fig. |^ for D = 3. The data collapse is obtained 
with the theoretical value 9 = 4/3. The scaling func- 
tion is composed of two distinct parts: a polydispersed 
small droplets distribution, with a small argument diver- 
gence of the scaling function associated with an expo- 
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nent r bigger than 1, well separated from a bell-shaped 
mono-dispersed distribution of bigger droplets centered 
around s = S(t). Most of the droplets in the system 
contribute to the small droplets distribution, which de- 
termines, since r > 1, the behavior of n(t), whereas the 
population of bigger droplets contains most of the mass, 
and S(t) is the typical mass of big droplets. The two dis- 
tinct populations of droplets are clearly visible on Fig. ||, 
which shows a typical configuration of droplets obtained 
by simulation of the deposition and coalescence model 
for d — 2 and D — 3 (spherical droplets on a plane). As 
shown in |sj the obtained droplets patterns are qualita- 
tively very close to the one obtained in some experiments 
of vapor deposition of thin films. 



power law behavior is well recovered in numerical simula- 
tion (see Fig. 0), and gives the results, r = 1.264± 0.002 




FIG. 3. A typical configuration in the scaling regime of de- 
position and coalescence of three-dimensional droplets on a 
plane. The picture shows a configuration after that 655360 
droplets of radius 0.75 were deposited on a 256 x 256 square 
surface with periodic boundary conditions, and S(t) has 
reached the value 10348.9. Two distinct populations of 
droplets are clearly visible: a monodispersed population of 
large droplets and a population of small droplets with a broad 
distribution. 



The exponent r can be determined from the numerical 
determination of the scaling function, but with important 
uncertainty due to statistical limitations. Thus, a better 
method is to extract r from n(t) oc (S(t)Y T ~ e \ This 



(D 



- 2), 
3). 



r = 1.18 ± 0.03 (D = 3), r = 1.074 ± 0.001 
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FIG. 4. The scaling law n(t) oc S(t) T ~ e for r > 1 is very 
well obeyed in numerical simulations in d = 1. The figure is 
a plot of log(n) versus log(S') for D = 3, for which we find 
t = 1.18 ±0.03. 



Direct determination from the scaling function yields 
r w 1.05 for D = 3, which is only in fair agreement with 
the previous result, but this underestimated value of r 
would be improved by using larger, or more numerous 
samples. The exponent r is strictly less than 9, but is 
nontrivial and cannot be determined by simple physical 
arguments. 

Therefore, it is interesting to derive a Smoluchowski 
equation for this model, and check the mean field value 
of t, if possible. Family and Meakin Q showed from scal- 
ing arguments that the coalescence kernel should have a 
homogeneity A = 2d/D — 1, but they did not determine 
its specific form. We proceed now to the determination of 
the equation. Neglecting multiple collisions, we examine 
the different events affecting the distribution N(s,t). 



Between t and t + 1, a droplet of radius s = kso is created as an outcome of the following processes: 

1. a droplet of mass sq falls on a droplet of mass si < s — sq, which occurs with probability Qi(sY D + s^ D ) d , Qi 
being a mass independent geometric factor. The droplet of mass si consequently reaches a mass si + so Then 
it coalesces with a neighboring droplet of mass s 2 — s — s± — so provided that they interpenetrate. 

The number of such events is, 



r!!ft 2 iV( Sl )jV( S2 )(sf +s,f) c 



(si+so) D -s 



G(s 1 ,s 2 ,r,t)(s{ > +s 2 D + r) 



d-l 



dr 



(2.5) 



(f^2 is another geometric factor). G(si, S2, r, t) is the probability density that a given droplet of mass si has a 
droplet of mass s 2 at distance s° + s£ + r as first neighbor. 
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2. a droplet of mass sq falls on a droplet of mass s — sq with which it coalesces, and the obtained droplet does not 
overlap any other droplet, 



number of events = CliN(s — sq)((s — sq) d + Sq ] 



s®-(a-a ) 15 



G{s - s , si,r,t)((s - s )o + s f + r) 



dr 



(2.6) 



3. a droplet falls in an empty space between droplets, 

number of events oc (1 — <j>(t))8 s>So , (2.7) 

where (1 — <p(t)) is the empty area fraction. 
A droplet of radius s disappears due to the following events: 
1. it coalesces with a droplet of radius s\ + sq which has grown 

1 A. 

number of events = 0irj 2 A r (s)A r (si)(sf o +s,f) d / G(si, s, r, t)(sf + s~ + r) d_1 dr. (2.8) 

Jo 



2. it grows: 



number of events oc N(s)(s D + Sq 



(2.9) 



To describe the large time scaling regime, we can take the continuous limit of small (but finite) so, to obtain the 
following continuous kinetic equation, 



d 1 f S 

d t N(s, t) + d s (s^N(s, t)) = - N{s u t)N(s - s 1 ,t)K(s 1 , s - S!,t)dsi 

2 Jo 



-N(s,t) / N(si,t)K(s,si,t)ds 1 +I(t)S(s-s ), 
Jo 



(2.10) 
(2.11) 



where the symmetric kernel K(x, y, t) is, 



K{x,y,t) = lim — / 
e Jo 



A Ax+e)d- x d 



G{x, y, r, i)(x D + y D +r) dr + symmetric 



(2.12) 



The time and mass units were redefined to eliminate multiplicative constants in the equation. I(t) is conse- 
quently renormalized to I(t) = c(l — <$>{t)), where c is a constant which could be easily determined, but is 
not essential to our discussion. It should be noticed that the distribution function is zero below sq at any t. 



I 

The mean-field approximation consists in neglecting is self-consistent: being proportional to the free surface 



spatial correlations i.e. in taking G(x,y,r,t) = 1. We fraction, it is a functional of N(s,t) (see Sec. HID) 
get: 



K (x, y, t) = (x D 1 + y D 'J^+K )" 



(2.13) 



This kernel has the homogeneity A = 2d/D — l derived 
from scaling arguments by Family and Meakin Eq. 
(2.10) is not a standard Smoluchowski equation since it 
incorporates two additional terms: an exogenous growth 
term d s (s iy N(s,t)), describing intercollision growth of 
droplets through absorption of small droplets, and a time 
dependent injection term. Moreover, the injection term 



I(t) is vanishing at large time, since the surface frac- 
tion covered by droplets goes to one in homogeneous 
growth processes, through renucleation in empty spaces 
(in heterogeneous nuclcation models, the coverage goes 
to a value <f> < 1). In fact, our numerical simulations in 
ID show that in the scaling regime 1 — (pit) oc n(t), as 
illustrated in Fig. [5]. This result is easily recovered from 
the scaling theory, 
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/+oo 
: 



0_ 

.-V) 



xnf(x)dx, (2.14) 



where Q is a geometric constant factor, which im- 
plies, since <p(t) — > 1 and = 1 + d/D, that 
J °° x d ' D f(x) = Or 1 . Since /(a;) oc x~ r , we see that, 

(1-0CO) oc (ao/S(t)) 1+,i/I, ~ T , if r > 1, and (l-0(t)) oc 
\s Q /S{t)) d / D , if r < 1, which yields (1 - ${£)) oc n(t). 
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FIG. 5. Plot of the free substrate area 1 — cj>(t) versus 
the number of droplets n(t) for a simulation of homogeneous 
growth with d — 1 and D — 3. At large time (small n), the 
two quantities are proportional, as understood in the frame- 
work of the scaling theory. 



B. Heterogeneous nucleation 
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FIG. 6. A typical configuration in the scaling regime of het- 
erogeneous nucleation of droplets in d = 2. The simulation 
was performed for lo — — 1 on a 512 x 512 square surface with 
256 2 droplets of radius 0.7 in the initial condition. 



Heterogeneous nucleation [Q corresponds to the case, 
common for water vapor condensation, when impurities 
on the substrate play a major role in droplet nucleation. 
A daily life example would be water condensation on a 
dusty pain. Nucleation occurs only on some nucleation 
centers, existing droplets grow from vapor, and coalesce 
when coming into contact, but no new droplet can nucle- 
ate in empty spaces. 

In the growth and coalescence model introduced by 
Family and Meakin one starts from an initial pop- 
ulation of droplets of same radius with no overlap. In 
the dynamics, individual droplets grow between collisions 
with, 



r = Ar u , 
or, equivalently (s = r D ), 

s = OAs' 3 , 



(2.15) 
(2.16) 



with j3 = (lo + D — 1)/D. In the following theoretical 
discussion we shall always set A = 1/D, but in numeri- 
cal simulations A was set to 1 (this just corresponds to 
a change in the time unit). We must have lo < 1 (or 
equivalently /3 < 1), otherwise the system is gelling as 
the mass of an individual droplet growing without colli- 
sion diverges at finite time. One step of simulation con- 
sists in increasing the radii of all droplets according to 



a discretization of Eq. (2.15), then tracking down and 
resolving all the resulting coalescence events, with the 
same rules as for homogeneous nucleation. 

Fig. |^ is a snapshot of a typical simulation of growth 
and coalescence of three dimensional droplets (D = 3) 
on a plane (d = 2). The simulations were carried out 
for lo — —1 on a 512 x 512 square surface with periodic 
boundary conditions from an initial population of 256 2 
droplets of radius 0.7 randomly placed without overlap. 
Qualitatively similar results are obtained when varying 
lo Pj. The droplets configuration in Fig. |^ appears to 
be visually very different from the one in Fig. ^ for ho- 
mogeneous nucleation. In heterogeneous growth, there is 
a single, monodispersed population of droplets. Another 
interesting feature of heterogeneous growth is that the 
surface coverage tends to a limit < 1 ( (|^J^] and ref- 
erences in B), which depends on D but not significantly 
on lo. The value of the asymptotic coverage was com- 
puted by Derrida et al. |Tc|| for several simplified models 
of coalescence in one dimension. Vincent f28|j , derived 
<f> = 0.57 from an approximate log-normal scaling solu- 
tion to a Smoluchowski mean-field equation (read below) 
with four-body collisions included for d — 2, D — 3 and 
lo = —2, in excellent agreement with the numerical value 
4> = 0.55. 

Consistently with Fig. ^, the droplets mass distribu- 
tion is bell-shaped H, as shown in Fig. [?] for lo = — 1. 
As in deposition and coalescence, an asymptotic scal- 
ing regime is reached at large time, as assessed by the 
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data collapse obtained for the mass distribution N (s, t ) 
at three different times. The scaling form of Eq. (2.1) 
was used with the theoretical value 9 = 1 + d/D, de- 
rived, as in the homogeneous case, from the fact that the 
surface coverage tends to a constant. 
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FIG. 7. A typical scaling function for heterogeneous nu- 
cleation with d < D. The figure shows the data collapse 
of N(s,t) in a simulation of growth and coalescence with 
oj = — 1, d = 2 and D = 3. The scaling function is clearly 
monodispersed and vanishes at xo ~ 0.2. 



The exponents z can be determined from the fact that 
S(t) is the only mass scale in the asymptotic regime, 
which implies that the distance between droplets n. 1 '"* 
scales as S 1 ^ . Then, from a rough evaluation of the to- 
tal collision rate, it is justified in || that the growth law 
of the typical droplet mass in the asymptotic regime is 
the same as that of an individual droplet in the absence 
of collision, except for a multiplicative constant renor- 
malizing the growth rate, 



which leads to, 



D 



\-(5 



(2.17) 



(2.18) 



These scaling results will be established for the corre- 



sponding Smoluchowski equation in Sec. Ill 



A corollary of this scaling behavior is that the scaling 
function f(x) cannot be diverging at small x. A simple 
argument even shows that f{x) is strictly zero below a 
finite x a > 0, as can be seen on Fig. |^ (see also Fig. 
2 in 1 10 1). In homogeneous nucleation, the lower cut- 
off of the distribution is the mass sq of the deposited 
droplets, which is time-independent, while S(t) diverges. 
Consequently, the lower cut-off of the scaling function is 



= \imso(t)/S(t), and f(x) can have a power law di- 
vergence when x — > 0. The situation is very different in 
heterogeneous nucleation. Consider the smallest droplets 
surviving at t. These are descendants of the droplets in 
the initial condition, that have not experienced any colli- 
sion since t — 0. As a consequence, the mass of the small- 
est surviving droplets sa(t) is, for a strictly monodis- 
persed initial condition N(s,0) oc S(s — sq), 

s (t)oc(*o(0) 1H3 + (l-/J)t)^ s (2-19) 

and So(t)/S(t) approaches a constant value xq > 0, in- 
dependent on so, when t — > oo. Since N(s,t) — for 
s < so(t), we see that f(x) = for x < x . Fig- 
ure ^ shows the scaling function obtained from simu- 
lations for d = 2, D = 3 and uj = —1. Numerical 
results give S(t) ~ 16t 1 /( 1 ~P\ with 8 — 1/3, whereas 
a (t) - ((1 - 8) Dt) 1 ^ 1 -^ from Eq. ( ^19|) (D appears 
in the formula since A = 1 in the simulation), which leads 
to xq « 0.2. This value of xq is fully consistent with Fig. 

B 

However, if ft = 1, because the growth of so(i) is ex- 
ponential in time, the situation can be different. The 
collisions renormalize the growth of S(t) and we expect 
at large time, 



S~(l + e(t))S(t), 



(2.20) 



where e is strictly positive, and might tend to zero, but 
there is no reason to expect that it should vanish so fast 
that S(t) oc so(t), and S(t) can be much bigger than 
So(i), leading to Xo = 0, and a possibly polydispersed dis- 
tribution. This behavior will be found in Sec. [II for the 
corresponding Smoluchowski equation. We shall see that 
S(t) is also much bigger than so(t) in mean-field in the 
case d = D, as S(t) oc (t lnf) 1 /( 1_/3 ', and that polydisper- 
sity can occur in this case. Actually, Family and Meakin 
H observed a qualitatively different, broader mass dis- 
tribution, with polydispersity, in numerical simulations 
for d = D. This case will be further investigated in Sec. 
IV. 

Now we proceed in deriving Smoluchowski's equation 
for this problem. In an early numerical and theoreti- 
cal work, Vincent [B8[ , considered heterogeneous growth 
of three-dimensional droplets on a two-dimensional sub- 
strate. His numerical simulations concerned early stages 
of growth (because he could not reach the large time 
asymptotic regime), when the number of droplets de- 
cays exponentially in time, in contrast with the power 
law behavior in the scaling regime, but he wrote down 
a mean-field Smoluchowski equation for heterogeneous 
growth. However, we shall find that his equation is in- 
correct since it does not conserve the number of droplets 
when the collision term is suppressed. 

Indeed, if we assume that droplets do not coalesce, we 
can find the contribution to Smoluchowski's equation due 
to the growth of individual droplets. The correspond- 
ing term must conserve the number of particle, since no 
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new droplet is introduced in the system: actually, the 
equation is just a continuity equation for the distribu- 
tion function N(s,t), and we find, 



d t N(s,t) + d s (s^N)(s,t) = 0. 



(2.21) 



If we bring coalescence into the picture, we note, fol- 
lowing Vincent pgfl , that the rate of coalescence of two 
droplets of masses si and S2 is, under the mean-field 
assumption, the time derivative of the cross-section oc 
(s\ + S2^ D ) d , which is proportional to (jj 1 dsi/dt + 
s y D ~ 1 ds 2 ldt)(s 1 ^ D + s^ )^ 1 . So that eventually, 



d t N(s,t) + d s (s l3 N)( S ,t) = 
\ f° N(s u t)N(s - s u t)K(si,s - ajdat 
N(s 1 ,t)K(s,si)ds 1 , 



N( S ,t)f + ~ 



with, 



K(x,y) = (xo +yo)(xo +yo) 



d-l 



(2.22) 



(2.23) 



Once again, redefinition of the time and mass unit was 
used to set multiplicative constants to one in the equa- 
tion. 

Vincent pq ], considered the case d = 2, D = 3, 
u> = — 2, relevant to the growth of epitaxial films. He de- 
rived Smoluchowski's equation for the radius distribution 
ip{r^t) = Dr 1 ^ 1 / N(r D , i), and found the correct colli- 
sion kernel, but his equation does not conserve the num- 
ber of particles when the kernel is set to zero, i.e. when 
collisions are ignored. The reason is that he erroneously 
derived that the change in ip due to growth alone was 
—r~ 2 d r ip, instead of the correct —d r (r~ 2 ip). As a conse- 
quence, some additional unphysical droplets arc created 
by his equation. This might be one of the reasons why 
Vincent had to include three and four-body coalescence 
events in his Smoluchowski equation to recover correct 
values for the fraction of area covered by the droplets, 
but the incorrect right-hand side in his Smoluchowski 
equation may as well have only minor consequences in 
his approximate computation. 

An interesting case arises when lu = 1 + d — D, since, 
as noticed by Family and Meakin this corresponds 
to the growth rate of large droplets due to absorption 
of deposited small droplets in homogeneous growth. In 
this case, Smoluchowski's equation describing homoge- 
neous growth differs from the one describing heteroge- 
neous growth only by the injection term. Exponents 
9 and z are the same for both models, and numerical 
simulations |8||| show that the scaling function for large 
droplets in homogeneous growth is very similar to the 
whole scaling function of heterogeneous growth. 

Thus, both growth and coalescence, and deposition 
and coalescence, are described in mean-field by a gen- 
eralized Smoluchowski equation, with additional terms 



accounting for intercollision exogenous growth of parti- 
cles (droplets). Therefore, it is interesting to perform a 
general study of this equation (with a generic kernel), 
and to see if its scaling behavior is consistent with the 
numerical results for droplets nucleation. 



III. SMOLUCHOWSKI'S EQUATION WITH 
GROWTH AND INJECTION 

We consider the following generalized Smoluchowski 
equation, 

d t N(s,t) + d s (s^N(s,t)) = 
\ J N(s 1 ,t)N(s - s 1 ,t)K{s 1 ,s - s-l) dsx 
- N{s, t) J N(si,t)K{s, si) dsx + I(t)5(s - 1), (3.1) 

K{x, y) being a general homogeneous kernel with expo- 
nents A and /j, defined as in Eq. (|1.2|). 

The equation describes a set of particles or clusters 
which collide with a mass-dependent collision rate K, 
and grow between collisions with (see below), 

s = s 13 . (3.2) 

Besides, some small particles (monomers) are injected 
with the injection rate I(t), with the possibility that / 
be a functional of N(s, t), as found for deposition and co- 
alescence in Sec. ||. A discrete version of this equation 
without the monomer injection term and with a constant 
collision kernel has been investigated for < (3 < 1 by 
Krapivsky and Redner |2£j . We shall see below that their 
results are independently recovered as special cases of our 
general discussion of the continuous equation, in the scal- 
ing regime where the discrete structure of the equation 
plays no role. 

When the growth and injection terms are absent, the 
equation reduces to Smoluchowski's equation, and its 
scaling properties have been extensively studied , 
but it is in no way trivial. Even in this case, very few ana- 
lytical solutions of Smoluchowski's equation are available. 
For the constant kernel K (x, y) = 1, an exact solution is 
known with N(s,t) - 4/i 2 e~ 2s /*. Other solutions 
concern the kernels x + y and xy |5l| . Despite its ap- 
parently simple structure, Smoluchowski's equation is yet 
another example of a highly nontrivial mean-field theory. 

In the following, we shall study the large time proper- 



ties of the solutions of Eq. (3.1), and we shall exhibit a 



rich diversity of behaviors depending on the parameters 
(3 (characteristic of exogenous growth) , A and fi (charac- 
teristic of the collision kernel). 



A. Gelation criterion 

A first interesting question is the possible occurrence 
of a gelation transition for such equations. Gelation cor- 
responds to the formation of an infinite cluster at finite 
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time. Without a growth term, nongelling kernels corre- 
spond to A < 1 |24|,^5|]. How is this modified ? In the 
absence of an infinite cluster, the evolution equation for 
the total mass in the system Mi (t) is obtained by multi- 
plying Smoluchowski equation by s and integrating over 
all masses, 

M 1 (t) = M p (t)+I{t), (3.3) 

which is physically obvious from s — s 13 . To discuss gela- 
tion, we have to be more cautious. Adapting the argu- 
ment for standard Smoluchowski's equation |3^j33|,|24|,|25|] , 
let us consider the mass transfer from clusters of masses 
s < L towards clusters of masses s > L, 

J L (t) = - sd t N(s,t)ds + aN(s,t)ds + I. (3.4) 
Jo Jo 



From Eq. (|3.l| ), we get 



J L (t)=L 1+ ^N(L,t) 

»L />+oo 

dxxN(x,t) / dyK(x,y)N(y,t), (3.5) 



J L — x 

where the first term is the mass flux through s = L 
due to the growth of individual particles, while the sec- 
ond term is the mass flux due to collisions. If there is 
no gelation Jl(£) must vanish when L — ► oo, and Eq. 
( |3.3| ) holds at any time. If there is gelation at t = t g , 
there is an infinite cluster, or gel, in the system for 
t > t g , and Jh{t) is nonvanishing for t > t g . At the 
gel point, Joc(t) = liniL^oo Jr,{t) may be infinite, but 
not for t > t g . The post-gel distribution must have a 
slowly decaying large s tail in order that Joo(t) be finite. 
If we make the ansatz N(s, t>t g )~ A(t) s~ T at large s, 



J 



L 1+f3 N{L,t) ~ A{t)L 1+l3 - T 
dxxN(x,t) / dyK(x,y)N(y,t)~A(t) 2 L 3+x - 2T dx dy xK(x,y)(xy)- T . 

ii J L-x JO Jl-x 



(3.6) 
(3.7) 



r 



We see that if gelation occurs, r must be equal to 
max(l + (3,(3 + A)/2). In the post-gel regime, the total 
mass contained in the sol phase (i.e. finite mass clusters) 
must be finite, which imposes r > 2. We conclude that 
no gelation occurs for max(A,/3) < 1. Anyway, (3 > 1 is 
forbidden since it leads to explosive growth of individual 
particles. If j3 = 1, Eq. (f]| yields M x {t) = e*, the total 
mass growth is faster than any power of t, which is the 
smoking gun of the gelling-nongelling boundary. 



This equation is easily solved for Z. With n(Q) — 1 and 
Z{z, 0) = Z (z), we find that Z(0, t) = n{t) = 2/(t + 2) 
and, 



Z{z,t) 



(* + 2 ) 2 (4^0(2) _ \ So (t'+2)i dt ') 



(3.10) 



which leads in the scaling regime t — ► 00, and for a 
monodispersed initial condition N(s,t) = S(s — 1), to, 



B. No injection 

For a while, we specialize to the case 1 = 0, corre- 
sponding to growth and coalescence. We first exhibit 
two exact solutions, then we make a complete study of 
the scaling solutions of the general equation. 



1. Exact solutions 



We can solve Eq. (3.1), in the case K{x,y) = 1, for 
(3 = and (3=1. To do this, we consider the Fourier- 
Laplace transform Z(z,t) of N(s,t), 



Z{z,t) 



+ OC 



5 N(s,t) ds, 



(3.8) 



Z(0,t) being the total density of clusters n(t) 



For (3 = 0, the Laplace transform of Eq. (3.1) reads 
1 



d t Z + zZ= ^Z 2 - Z(0,t)Z. 



(3.9) 



N(s,t) 



t 2 Int 

The total mass in the system is, 



(3.11) 



Mi(t) = -d z Z(z = 0,t) = l + 21n(i+ -). (3.12) 



For (3=1, the equation for Z is, 

d t Z - zd z Z = ^Z 2 - Z{0, t)Z, 



(3.13) 



and once again n(t) = Z(0,t) = 2/[t + 2). If we choose 
the variable u = ze , the equation reduces to a first order 
differential equation in time, and the solution is, 



Z(z,t) 



2Z (ze t ) 



t + 2(l-Z t 



(ze t ))(t + 2) + 2Z (ze t )- ^ 



With an initial monodispersed distribution, Zq(z) = e~ z , 
Z(z, t) has a pole at Zo(t) = ln(l + 2/t), and we can 
explicitly compute N(s,t), 
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N(s,t) 



{t + 2) 2 e> 



■ exp (-se _t ln(l+2/t)), (3.15) 



which leads in the large time limit to, 



t z e l 



(3.16) 



The corresponding solutions for the discrete Smolu- 
chowski equation have been independently derived by 
Krapivsky and Redner and coincide with the so- 

lutions above in the scaling limit. This coincidence 
was to be expected since, in the scaling regime, the di- 
vergence of S(t) leads to the oblivion of the discrete 
structure of the equation. This provides a solid con- 
firmation of the surprising result that, although the 
growth term did not change the scaling function, still 
equal to f(x) = 4e~ x , the scaling is no longer of the 
form N(s,t) ~ S(t)- e f(s/ S(t)). Rather, the scaling is 
N(s,t) ~ Y(t)- l f(s/S(t)), where Y(t) oc S 2 {t)/M 1 {t) is 
not a power of Sit). This form enforces n(t) oc l/t, a 
result that holds for any f3 between and 1. 

A consequence of the logarithmic correction in the case 
f3 = 0, is that in contrast to what we stated for the generic 
case in Sec. 



[IB 



Xq is equal to zero. The reason is that 
S(t) in this case grows faster than do individual particles 
in absence of collisions. Hence So(t)/S(t) oc 1/lni goes 
to zero. This point will be fully discussed in the case of 
a general kernel under the dynamic scaling assumpt ion. 
For (3 — 1, we see that the scenario discussed in Sec. [I B 
occurs, S(t) ~ te* corresponds to a slowly vanishing e(t) 
in Eq. ( |2.20|) . 



2. Scaling theory 

We would like to study the scaling properties of Smolu- 
chowski's equation. Even though Smoluchowski's equa- 
tion results from an approximation, its scaling behavior 
is usually exact above an upper critical dimension d c , and 
is in many cases qualitatively correct even below d c . 

Some simple arguments may give a qualitative under- 
standing of the different regimes to be expected for the 
equation. Indeed, if we suppress the collision term (i.e. 
the right hand side) in Eq. (3.1), we are left with a conti- 
nuity equation describing a set of particles which grow in 
time with s = s@, and is associated with the mass scale 

Sg{t) OC iV(l-/3). 

Conversely, if we suppress the exogenous growth term 
d s (s@N) in the left hand side, we are back with a stan- 
dard Smoluchowski equation describing clustering with 
mass conservation. The scaling properties of this equa- 
tion are well-known pj-pq]. The typical mass in the 
scaling regime is S c (i) oc F'( 1 ~ A ), and 6 — 2. 

Thus, when both exogenous growth and collisions are 
active, we expect to observe a "competition" between 
the two dynamic mass scales S c and S g . If f3 < A, 



S g (t) 3> S c (t), and in the scaling regime we expect 
S(t) oc S g (t) and z = 1/(1 — /?)■ If f3 > A, on the contrary, 
the typical mass of particles increases essentially due to 
collisions, hence S(t) oc t 1 /^ 1- *) an d z = — A). In 
the marginal case A = (3, logarithmic corrections to S(t) 
may be observed. In fact, we know from the exact solu- 
tion ofif=l,/3 = = A that such corrections actually 
occur. This leads us to a slightly more general scaling as- 
sumption than the one we made for droplets coalescence 
models, 



*(M)~y(t)-v(^). 



(3.17) 



We do not assume a priori that Y(t) oc S(t)~ e , since we 
know from exact solutions that it is not always true. 

Notice that the scaling function f(x) is not uniquely 
defined by Eq. (3.17) unless we give a precise definition 
of S(t) and Y(t). For instance, a widely used definition 
of S(t), (and the one actually used in numerical simula- 
tions), is 



S(t) 



< s > 



(3.18) 



However if we know a scaling function / s (a;) for given 
definitions of Y and S, any other scaling function f(x) 
corresponding to other definitions, is related to f s by 



f(x) = Kf s ((x), 



(3.19) 



k and £ being two constants. 

From the picture above, it is obvious that the physical 
cut-off, i.e. the mass So(t) below which N(s,t) is strictly 
zero, scales as S g (t). This is just the translation in terms 
of Smoluchowski's equation of the discussion we had for 
droplet growth and coalescence. Since S(t) > s a (t), ei- 
ther S(t) and S g (t) have the same scaling and the scaling 
function f(x) is zero below a certain argument x$ > 0, or 
S(t) ^> S g (t), xq is equal to zero, and / may have a small 
x divergence f(x) oc x~ T , with a polydispersity exponent 

T > 0. 

The scaling of the moments of the distribution N(s, t) 
is altered by the existence of a polydispersity exponent. 



s a N( Sl t)ds — I _ x a f{x)dx. 

so(t) 



Y 



■n't' 
-'it, 



(3.20) 



If there is no polydispersity exponent or if t < 1 + a, the 
integral tends to a finite limit when t — > oc, and, 



M a (t) oc 



S 1+a 
Y 



If t > 1 + a, the integral diverges and, 
M a (t) oc ^ So (t) 1+a - T . 



(3.21) 



(3.22) 
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Finally, if r = 1 + a, 



> , , N Sl + a , S (t) 

M a (t) cc— In ^. 



(3.23) 



Under the general scaling assumption, we get the fol- 
lowing scaling for the different terms of Smoluchowski's 
equation: 



Another important equation is Eq. (3.3) for the evolution 
of the total mass in the system, which, in the absence of 
injection, reduces to, 



Mi = Mp. 



(3.27) 



9tN(s,t) ~ -~ !yf(x) + ^xf'(x)\ , (3.24) 



qP-l 

d.(80N(8,t))~Z—(xPfy(x), (3.25) 



s 1+x 

Collision term - -rrzz- (...). (3.26) 

y A 



Now, it is possible to find the asymptotics of Mi, F 
and S*, depending on the values of A and 0, under the 
sole scaling assumption. In fact, though the line followed 
in the demonstration is quite simple, details are rather 
intricate due to the multiple cases to be examined. A 
full length discussion is given in Appendix and results 
are summarized in Table |. Here, we shall only comment 
some interesting points. 



parameters 



Mi(t) 



S(t) 



Y(t) 



1 > > A 



Mi cx S f) - X 



S(t) cx t z 
1 

z 



1-/3 



KocS* 
9 = 2 + X 



1 > A > 



Mi 



constant 



1 - A 



= 2 



1 > X = 



hit, if/x<0 
(lni)ln(lni), if fx > 



S(t) cx (tMi) 2 
1 



c2 

Mi 



= 1 
< A < 1 



Mi = Mi(0)e* 



S(t) cx e~ 



Y(xS 2 e 



= 1 
A = 



Mi = Mi(0)e* 



cx te* 



Y ottS 



A = 1 > 
fi>0 



Mi — > constant 



5(t) cx e foV * 



y cx s 2 Vt 



\ = l> 

M < 



Mi — > constant 



cx e fct 



Fcx^ 



A = = 1 
H < 



Mi = Mi(0)e* 



S(f) 



Y cx 5 2 e _t 



A = /3 = 1 
/i > 



Mi = Mi(0)e* 



S(t) cx e b ^ 



Y cx S^e"* 72 



TABLE I. Results of the scaling theory 



The scaling theory is consistent with the qualitative 
discussion above based on the idea of competing dynam- 
ical scales. It is found, that for A < < 1, S(t) scales 
as s {t) ~ t 1 /(i-/3) j Y (t) cx S{t)- , with 6 = 2 + A - 0, 
and the scaling function is zero below a finite xq. If wc 
come back to droplets models, this A < condition just 
corresponds to d < D, and we find 9 = 1 + d/D. Hence, 
the scaling results of the mean-field theory are in full 



agreement with the discussion and results in Sec. [I B 



For A > 0, S(t) scales as t 1 ^ 1 ^^, the mass is asymp- 
totically conserved with 8 = 2. The scaling function may 
have a polydispersity exponent, since now S(t) ^> so(i), 



and the scaling equation is, 

b [xf(x) + 2f(x)] = f(x) / +o ° f( Xl )K(x, Xl )d Xl 

_ \ Jo f( x i)f( x ~ xi)K{xi,x- x\) dxi, (3.28) 

that is the same scaling equation as for standard Smolu- 
chowski's equation without growth of particles, which 
makes it possible to use all the corresponding results or 
techniques [E4|-p6| . This case corresponds to d > D, but 



as further discussed in Sec. IV, heterogeneous growth is 
always gelling in this case, and the mean-field approxi- 
mation breaks down. 

For A = 0, we find that S{t) is no longer a pure power 
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law, but incorporates logarithmic corrections. The total 
mass in the system increases logarithmically, Mi (t) oc In t 
and S(t) ~ (tMi(t)) 1 /(i-«. Once again, S(t) > s (t) 
and the scaling function is Eq. (3.28). Thus, there is a 
polydispersity exponent if the kernel has fi > (see be- 
low), and there is an addition ln(lnt) correction for fi > 
kernels. For heterogeneous growth, A = (3 corresponds to 
d = D, and the mean- field theory accounts for the qual- 
itative difference between d = D and d < D observed in 
numerics (see 8| and below). This point will be fully dis- 
cussed in Sec. [V. This also recovers the scaling behavior 
of the exact solution for K = 1 and (3 = 0. 

For (3 = 1, the scaling of the exact solution K = 1, 
(3 = 1 is recove red. For A > 0, the scaling equation is 



once again Eq. ( 3.28 ). Other results in Table |, show the 
great diversity of scaling regimes depending on (3, A, and 

For a constant kernel, and 1 > (3 > 0, we recover 
the result of Krapivsky and Redner (2^] . The latter au- 
thors assumed that the scaling function has essentially 
the same shape as in the case (3 — or the pure ag- 
gregation case. From our analysis, we know that this 
assumption is actually not verified. However, it can be 
seen that the key point of their demonstration is that 
f(x) has no small x divergence, which is indeed true. 

The fact that the scaling results of the d < D growth 
and coalescence are recovered by Smoluchowski's equa- 
tion approach gives a firm basis to the heuristic argu- 
ments used to find a posteriori the exponents from the 
obtained numerics. Moreover the kinetic equation ap- 
proach is predictive, and provides a synthesized classi- 
fication of the aggregation models, depending on a lim- 
ited number of relevant parameters, provided that the 
approximation be justified. 



3. Polydispersity exponents 

An interesting corollary of the scaling theory of gen- 
eralized Smoluchowski's equation with growth, is that in 
the cases (3 < A and (3 = 1, the scaling equation is ex- 
actly the s ame as for standard Smoluchowski's equation, 
Eq. ( |3.28| ), where b is called the separation constant, 
which we set to one by absorbing it in the scaling func- 
tion (which corresponds to a redefinition of Y). Note that 



if / is a solution of Eq. (|3.2 



„1+A 



f(cs) is also a solu- 



tion, which corresponds to different possible definitions 
of S(t) (remember the discussion above Eq. (3.19)). 

Thus, all the results known for the scaling function of 
standard Smoluchowski's equation also hold for the gen- 
eralized equation. For instance the scaling function f(x) 
of the K — 1, (3 — 1 case can be derived from the exact 
result for standard Smoluchowski's equation with K = 1, 
for which fo(x) = e~ x is a scaling function. For a given 
definition of S and Y, the corresponding scaling func- 
tion for K = 1, (3 = 1 is obtained using f s = fo in Eq. 



( p^| ). If we use Eq. (|1J|) as a definition of S(t), £ is 
constrained to the value £ = 2 and if we define Y(t) by 
Mi(i) = S 2 /Y, we get n = 4, which leads to, 



;vM)°c§^ 2s/s , 



(3.29) 



with S(t) oc te l in agreement with the exact result Eq. 

(Hi). 

We can also find the exact scaling function for (3 = 1 
and K{x, y) = x + y (which corresponds to A = 1). For 
standard Smoluchowski's equation, a scaling function is 



-3/2 



e~ x . Thus, it is also a scaling function for 
Eq. (3T), and we obtain the e xact result r = 3/2. 

The scaling equation Eq. ( |3.28| ) for a general kernel 
was extensively studied in the literature. Van Dongen 
and Ernst [p4|]25| l showed that the qualitative shape of 
the scaling function f(x) at small x depends on two pa- 
rameters, the homogeneity degree of the kernel A and the 
exponent /j defined by Eq. ( |l.2j ). 

For n < 0, the scaling function vanishes as exp(— ax fl + 
o(x fJ ')) at small x, and there is no polydispersity expo- 
nent. 

For kernels with /i > 0, there is a polydispersity expo- 
nent r = 1 + A. For /i = 0, there is also polydispersity, 
but with a nontrivial exponent r < 1 + A, 



x x f(x)dx. 



(3.30) 



The determination of t for /i = was a challenge un- 
til recently, because solving numerically Smoluchowski's 
equation proved rather heavy, and often unsuccessful. 
For the most studied [i = kernel, 



Ki(x,y) = (x° +yvy 



(3.31) 



which corresponds for instance to Brownian coalescence 
(in a d + 2 dimensional space) with mass-independent 
diffusion constant, very few values of the exponent r 
were known. In d — 1, Kang et al. (34), noticed that 
direct numerical resolution of Smoluchowski's equation 
did not reach the actual scaling regime, but a pseudo- 
asymptotic regime, with apparent scaling but wrong ex- 
ponents (some exact bound were known for r ]35| , p6| ] ) . 
For D — 1 and d < 1, the scaling regime could be reached 
by Krivitsky J36|, leading to the determination of r for 
ten values of d between and 1. Song and Poland [ B7| |, 
used a power series in time expansion method, to treat 
the cases (d — 1,D = 2), and (d = 2, D = 3), but we 
showed |2(| that their result in the latter case is in con- 
tradiction with some exact bounds. 

We recently introduced a general variational method 
for computing accurate values of r at very low numeri- 
cal cost p6j. This method was used to make a complete 
study of the polydispersity exponent of the kernel, 
for a wide range of d and D. We chose to start directly 
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from the scaling (infinite time limit) equation Eq. (3.28). 
We did not try to solve this equation for the whole scal- 
ing function, but focused on the determination of r. The 
key relations which we used are, on the one hand, inte- 
gral equation Eq. (3.30) for r, and on the other hand a 



series of integral equations for the moments M a of r |24| ] 
obtained by multiplying Eq. ( 3.28| ) by x a and integrating 
over x, for any value of a > r — 1 (such that the integrals 
converge in zero), 



2(1 -a) / x a f{x)dx 
Jo 



f(x)f(y)K(x, y) [x a + y a - (x + y) a ] dxdy. 



(3.32) 



Combining these equations one obtains, 



r = 2- (1 - a) 



//o°° g( x ,v) dxd v 

ffo°° 9(x, y)A{x/ y) dxdy ' 



(3.33) 



with g(x,y) = f(x)f(y)(x a y x + x y a ) and A(u) = 
(1 + u a - (1 + u) a )K(l, u)/(u a + u x ). 

The ratio in Eq. ( |3.33| ) is the inverse of the average of 
A(x/y) with weight g(x 7 y), so that computing the maxi- 
mum M a and the minimum M a of A for various values of 
a leads to exact bounds for r that can be used to check 
nume rical evaluations of r (sec ||(| for details), since Eq. 
( 3.33 ) implies, 



(1 - a)/m a < t < 2 - (1 - a)/M a . 



(3.34) 



The idea of the variational approximation is to choose a 
parametered family of variational functions, and to min- 
imize the violation of Eq. ( [3.33 ) for a well-chosen sample 
of values of a. The key point is the choice of the varia- 
tional function. As argumented in p^| , a natural three 
parameters class of functions, is: 



fv(x,T ,Ci,c 2 ) 



1 

x r ° 



Cl 



C2 



£1-1(7-0) x > 



(3.35) 



The last term corresponds to the exact asymptotic decay 
at large x of the scaling function |2jj2^], while To is the 
polydispersity exponent, and t\ is the subleading expo- 
nent in small x (its value in function of To is taken to 
be the same as for the exact scaling function, see below 
and (26|). This class of function has the correct large x 
and small x asymptotic behavior expected for the scaling 
function. Besides, it contains the exact scaling functions 
for K — 1 and K = x + y, therefore the variational ap- 
proximation yields the exact result for r in these cases 
(as checked in p6fl). 

W e use d as error function, measuring the violation of 
Eq. ( 3.33 ) for a set of n moments on , 



X 2 (fv) = ^2(ro-G ai (f v )Y 



(3.36) 



where G ai (f v ) is the right hand side Eq. (3.33) for 
a = oii and f = f v . This error function is by construction 
strictly zero for the exact scaling function. For the chosen 
class of variational functions, G ai (f v ) can be expressed 



in terms of Gamma functions and simple II? integrals, 
which makes its numerical computation extremely fast. 

In p^j, the variational method was used to perform 
a complete study of the kernel Kp, for < d < 3 and 
d < D < 7. Some analytical expansions were obtained for 
r in the limit d — > , D — > 00 and d — > 00 with A = dj D 
constant. These analytical results combined with exact 



inequalities obtained from Eq. (3.33), were used to check 
the variational results. The obtained results were in ex- 
cellent agreement with the few existing numerical values 
for r as well as with the asymptotic expansions. Due to 
its extremely low computational cost of the variational 
approximation, compared to other methods in the lit- 
erature, and its excellent accuracy, the variational ap- 
proximation seems to be a good practical solution to the 
problem of the determination of t. 

From this study of generalized Smoluchowski's equa- 
tion with growth of particles, we see that the variational 
method can also be used to find the r exponent for this 
equation, when polydispersity occurs, i.e. when A = ft 
or f3 = 1. An interesting physical application of these 
results is heterogeneous growth with d = D, for which 
Smoluchowski's equation derived in Sec. II B is in the 
class P = A. 

For this problem, (3 is equal to 1 + (u> — 1)/D and the 
kernel is, 



K(x, y) = (x D +y D )(x D +y D ) 



(3.37) 



This kernel is formally similar to the one describing diffu- 
sion limited cluster-cluster aggregation p8|-^0|, but the 
meaning of the parameters is different. For this kernel 
we have, 



0, if w > 
uj/D, if lo < 



(3.38) 



From the scaling theory we expect that S(t) oc (t\nt) z 
with z = D/(l — bS), and a transition from a polydis- 
persed scaling function with a nontrivial r exponent, for 
lo > 0, to a small x vanishing scaling function, for lo < 0. 

Consequently, it is interesting to determine the mean- 
field polydispersity exponent r for this kernel using the 
variational approximation. This will be done in Sec. |rvj , 
in which we also present comparisons of scaling results 
from Smoluchowski's equation approach with direct nu- 
merical simulations of heterogeneous growth with d = D. 
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C. Constant injection 



z(Mp + I + z^ 1 - 13 ) oc (z 2r - 2 - A ), (3.44) 



We now turn to the case of a constant injection rate. 
Interest in aggregation models with injection was origi- 
nally aroused from applications in chemical engineering 
(coagulation in stirred tank reactors) and atmosphere sci- 
ences |l^,^p|JT4]|4^Jl^8| . In these contexts, injection 



was often associated with a sink term. The emergence of 
the concept of self-organized criticality [ ^3|j44| resulted in 
a renewal of interest in aggregation models with constant 
injection |L9]-[2l]], since these systems commonly evolve 
to a steady state asymptotic power law distribution, and 
therefore provide examples of self-organized critical sys- 
tems. This behavior is assessed by numerical simulations 
and exact solutions in one dimension ]2(],^lJ . 

Hayakawa p^ j studied Smoluchowski's equation with 
injection of monomer. He showed that for non gelling sys- 
tems, with A < 1 p4| , the asymptotic steady state had a 
power law large s decay with an exponent r = (3 + A) /2. 
Here we shall investigate the steady state in the presence 
of a growth term with exponent 0. We shall suppose for 
convenience that the coagulation kernel K(x, y) is equal 
to x^y v + x v y^ . A = fi + v is the homogeneity degree of 
the kernel. The results are however true for any kernel. 

We are interested in the asymptotic steady state 
reached by the system at large time. We shall see that it 
has a large s power-law decay with an exponent r that 
we are able to compute in terms of A and 0. To achieve 
this program, let us call Z a (z,t) the Laplace transform 
of s a N(s, t) defined by: 



p+oo 

Z a (z,t) = / s a N(s,t)e- zs ds, 
Jo 



for which we get the following equation: 



(3.39) 



d t Z 1 + zZ & = Z^Z, y - Z„M V - ZvM^ + Ie~'. (3.40) 

Now we consider the equation for the steady state, 

zZf = (Z™ - M~)(Z~ - M?) + I{e~' - 1). (3.41) 

The large s behavior of the steady state distribution is 
reflected in the small z behavior in Laplace space, 



ZZ>(z)-MZ>~ Ca z T <*, 



(3.42) 



if M 1 ^ is finite. If is infinite, it does not appear in 
the left hand side. Note that is certainly infinite 

because there is constant injection of monomers and no 
dissipation of mass (at finite time). As a consequence 
< r a < 1 for all a < 1. If r Q is non integer then for 
s — > oo, 



s a N(s,t = oo) ~ -^L(l + r Q )s 1 - T °. 
As a consequence, t — 1 = tq = r a + a, and, 



(3.43) 



hence if r-1-0 > 0, then 1 = 2r-2-Ai.e. r = (3+A)/2, 
whereas if r — 1 < 0, Mg° is infinite and does not appear 
in the left hand side of Eq. fl3.44| ) then t-/3 = 2t-2-A 
i.e. t = 2 — + A. To summarize, we find, 



(3 + A)/2, if 0< (l + A)/2 
2 + X-0, if 0> (l + A)/2 



(3.45) 



Thus, we see that the exogenous growth term brings 
in a new feature: above a critical growth parameter 
f3 c = (l + A)/2, the power law exponent of the asymptotic 
state depends continuously on 0, whereas if (3 is less than 
C , the exponent is unaffected by the growth term. 

The case (3 — (1 + A)/2 requires some additional care. 
Interpolation of the two regimes above would lead to 
t = 1 + j3, and it is possible to show that there is a 
logarithmic correction N(s, +oo) oc l/(s 1+/3 lns). 

For (3 < 1 + A, t has the value found by Hayakawa 
p8[ in the absence of exogenous growth from the same 
Laplace transform arguments. However, we find that 
these demonstrations in Laplace space are not very il- 
luminating, as they do not really show what happens 
"physically". Here we would like to point out that this 
exponent can be directly found from a simple argument 
for N(s,t). For convenience let us first forget the exoge- 
nous growth term (the argument is the same), and let 
us consider the steady state condition. N 00 (s) is a sta- 
tionary distribution, in the sense that if we start from 
N(s,t = 0) = N OQ (s), then the distribution does not 
evolve. From this point of view, it becomes clear that 
the total mass injection rate I must exactly be compen- 
sated by the mass dissipation by collisions. Thus, the 
total mass flux due to collisions must be finite. This just 
means that the steady state is at a gel point, and the 
argument of Sec. Ill A can be readily adapted to obtain 
t = (3 + A)/2. Furthermore, as the total mass is infinite 
in the steady state, we do not have the restriction r > 2, 
which determines the gel criterion for gelation at finite 
time. Here, the transition occurs at infinite time, when 
the system has self-organized to the critical point of a gel 
transition. 

If we bring the exogenous growth term into the picture, 
we can also find the exponents from the same argument. 
Now, the mass injection rate is I + Ml 







If Mf is finite, 



we still find r = (3 + A)/2. If > (1 + A)/2, is 
diverging with such a value of r. Consequently, Mp is 
infinite, and the steady state condition is now, 



lim 



s /3 A^ 00 (s)ds-C(L) =0, (3.46) 



where C'(L) is the integral in Eq. (3.5). The vanishing 
of the divergence imposes 1 + 0— r = 3 + A — 2r, and we 
recover t = 2 + A — 0. 
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D. Constant mass injection. 

Now, we would like to turn back to the initial prob- 
lem of homogeneous nucleation, and the corresponding 
Smoluchowski equation. Let us remark that the colli- 
sion kernel has A = 2d/D — 1 = 2/3—1, just on the 
border line of the two re gime s for constant injection. 
However, we saw in Sec. II B that the injection term 
of small droplets vanishes as 1 — <j>. In the derivation 
of Smoluchowski's equation for droplets deposition and 
coalescence, we found I(t) from a geometrical argument. 
Another way of seeing it, directly from the Smoluchowski 
equation, is to impose the additional constraint to Eq. 
(3.1) that the mass injection rate be a constant (as in ho- 
mogeneous nucleation), say Mi = 1. From M\ = Mp + I, 
this is equivalent to I(t) = 1 — Mp(t). For droplet depo- 
sition, and with this choice of constants, Mp(t) = <f>(t), 
and the geometrical argument is recovered. 

Thus, we shall now discuss the case of A = 2/3 — 1 and 
constant mass injection Mi = 1, i.e., 



I(t) = l-Mp(t), 



(3.47) 



for (3 < 1. Once again, we shall make the scaling as- 
sumption of Eq. (3.17). As in homogeneous nucleation, 
Mi(t) oc t leads to Y(t) oc S 2 jt. A very interesting result 
is that Mp must tend to 1 at large time. First, it is easily 
seen that Mp(t) cannot diverge. The reason, is that the 
injection rate of "area" into existing particles is equal to 
(3M 2 p-i and is always dominated by Mp{t) in the scaling 
regime. More precisely, the evolution equation for t he oc- 
cupied area fraction Mp is obtained from Eq. (3T), and, 
since collisions cannot increase Mp ((3 < 1), we have the 
inequality, 



Mp<f3M 2p _i + 1-Mp(t). 



(3.48) 



Then, since (3 < 1 implies 2/3 — 1 < (3, it is possible to 
show that for any value of a possible polydispersity ex- 
ponent, we have in the scaling regime, [Mp — [3M 2 p-i) ~ 
cMp, where c is a strictly positive constant. In fact, 
Mp > M 2 ^_i and c = 1, if r < 1 + /3 while Mp oc M 2 p-i 
and c = l/(r -1-/3)- /3/(t - 2/3), if r > 1 + (3. This 
result, combined with Eq. 



(3. 



leads to, 



Mp < 1 — cMp(t), c > 0, 



(3.49) 



J 



which shows that Mp cannot diverge. It is also clear that 
there is no way that Mp could become negative, since the 
smaller Mp, the larger I(t). Therefore, if we rule out any 
pathological oscillatory behavior, Mp tends to a constant 
</>, that may possibly be zero. Now, if </> ^ 1, then the in- 
jection is asymptotically constant, and, from Sec. Ill C| , 
there is a critical steady state, with N(s) ~ l/(s 1+l3 In s) 
at large s, and Mp diverges, which is contradictory (be- 
sides, if <fi > 1, the distribution is negative near s = 1). 
Thus <j) ~ 1, which is a nontrivial result, and was a not 
obviously fulfilled necessary condition for our mean-field 
approach to correctly describe droplets deposition and 
coalescence. 

Now, let us discuss the scaling properties of the equa- 
tion. Since the injection term is vanishing, we expect 
the scaling equation (which describes large clusters) to 
be the same as in the case without injection. However, 
the fact that the cut-off sq is constant, and therefore neg- 
ligible compared to S(t) selects a solution different from 
the one obtained without injection. 

To be more precise, we know that Y ~ S 2 /t and that 
Mp has a finite limit. If we assume that there is poly- 
dispersity with t > 1 + (3, those two conditions lead to 
S(t) ^> i 1 ^ 1- ' 3 ', and we find that the scaling equation is 
once again Eq. ( 3.2S ), which yields r<l + A<l + /3, 
in contradiction with our assumption. Thus, r < 1 + (3, 
and Mp cx iS'(£) ,3 ~ 1 , leading to, 



1 



1-/3' 



(3.50) 



which corresponds to the results previously obtained for 
droplets deposition and coalescence (with (3 = d/D). The 
scaling equation is Eq. ( |A22| ), with positive a and b. This 
equation is nonlinear and is likely to admit several classes 
of solutions. We have seen that when there is no injec- 
tion, a solution is selected which vanishes below a finite 
xq > 0. However, in presence of injection the scaling 
function has no lower cut-off (xq — 0), and we can have 
a polydispersity exponent. 

To investigate the small x behavior of /, we introduce 
the auxiliary function (p(x) = x^^ 1 f(x), which leads to 
a new scaling equation, 



2s 1 - / V(s) + s 2 -V0) - scp'(s) - ip(s) = <p(s) / ip(si)K(s, si)ds 



-hoc 



tp(si)if(s - si)K(si,s - si)dsi 



(3.51) 



where K(x, y) = x 1 ~Py 1 ~°K(x, y) (e is included to regu- 
larize the collision terms which are separately diverging in 
the e — > limit f24|]). We remark that the most diverging 



term for x — > in the left hand side is — sip'(s) — ip(s), so 
that, as far as the determination of the asymptotic behav- 
ior ip(s) oc s~ r is concerned, we can straightforwardly 
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generalize the results of van Dongen and Ernst p4| , p5| . 
The kernel K has the homogeneity A = 2 + A — 2/3, and 
fl = 1 + /J, — ft, and we find that the function vanishes at 
small x for p, < 0, that r' = 1 + A for jl > 0, and, r' is 
non trivial, with, 



r+oo 

t' = 1 + / tp(x)x x dx < 1 + A, 
Jo 



(3.52) 



for ft = 0. 

Therefore, for fi < /3 — 1, we have no polydispersity, 
while for > ft — 1, we have r = 2 + A — ft, and for 
fx = ft-l, 

r + co 

t = /3+ f(x)x x+1 - )} dx < 2 + A - ft. (3.53) 
Jo 

Now, for A = 2ft - 1, wc find that r = 1 + ft, if 
/i > ft — 1, while r is nontrivial and strictly less than 
1 + /3 for fi = ft — 1. Hence, for // < ft — 1, the scal- 
ing theory is consistent, while for /i > /3 — 1 there is a 
contradiction with r < 1 + /3. The latter cas e pre cisely 
corresponds to droplets deposition (see Eq. 2.13 ), and 
we cannot conclude. However a consistent scaling with 
a nontrivial polydispersity exponent could be obtained if 
we include pair correlations into the collision kernel and 
if the resulting kernel has /i = ft — 1. 

This scenario is supported by an early numerical work 
of Tanaka |45| . Tanaka solved a set of coupled differential 
equations describing growth and coalescence with renu- 
cleation for d — 2, D — 3. Dynamical pair correlations 
due to excluded volume were included in an approximate 
form. It seems quite clear that his set of equation be- 
comes equivalent to a Smoluchowski equation very simi- 
lar to ours in the large time limit, but with a collision ker- 
nel modified by correlations, and Tanaka found a bimodal 
droplets mass distribution with a nontrivial polydisper- 
sity exponent, in agreement with the results described 
in Sec. |j It would be interesting to try to determine 
the kernel from his equation, although this seems quite 
difficult. 



IV. 



HETEROGENEOUS GROWTH WITH 
POLYDISPERSITY 



In Sec. Ill, we found from a mean-field approach that 
the kinetics of heterogeneous growth with d = D (for 
instance, discs on a plane, or spheres in 3D), should be 
qualitatively different from its counterpart with d < D. 
From the scaling theory of generalized Smoluchowski's 
equation, we found that there should be a transition from 
a monodispersed scaling function for w < 0, to a poly- 
dispersed function with a nontrivial polydispersity expo- 
nent r for lu > 0. This mean-field result is actually very 
interesting, since it corroborates numerical simulations 



performed by Family and Meakin [|3j|5j], who found that 
polydispersity occurs for d = D = 2 and lu = 0.5. 

Thus, our Smoluchowski equation approach sheds 
some new light on heterogeneous growth with d — D, 
which was not much studied due to the fact that interest 
was primarily focused on d = 2, D = 3 relevant to breath 
figures, and also to the fact that numerical simulations 
are much more difficult in this case (see below). In this 
section, we first fully discuss what should be expected 
from the mean-field theory, and compute the polydisper- 
sity exponents for u> > 0. Then we present some numeri- 
cal simulations in d = 2 and discuss the relevance of the 
mean-field theory. 



A. Mean-field theory 



In Sec. [IB, it was found that the collision kernel cor- 

D was, 



responding to heterogeneous growth with d 
K(x,y) = + y%)(x~5 + y~ 



■.D-l 



(4.1) 



with X = 1 + (io~ 1)/D — ft. The corresponding general- 
ized Smoluchowski equation was found to be nongelling 
for u> < 1, and in the following we shall take u < 1. 

The reason why growth with d = D is different in 
mean-field from d < D, is rather subtle. As discussed in 
Sec. 



Ill 



the "competing" dynamical mass scales corre- 
sponding to exogenous growth and growth by collision, 
respectively S g (t) and S c (t) are of the same order at large 
times for d = D, which leads to a marginal enhancement 
of the growth of the typical mass S(t), and to logarithmic 
mass growth, 



S(t) oc (t In t) — , Mi (t) ~ hi t. 



(4.2) 



This implies that the cut-off x = \im s o (t) / S(t) in the 
scaling function is zero, in contrast to the d < D case for 
which xq > 0, and the scaling equation is the same as for 
Smoluchowski's equation without growth. 

For lu > 0, we have fj, = 0, and consequently there is 
a nontrivial polydispersity exponent r. We can use the 



methods discussed in Sec. Ill B 3 to study r 



As a preliminary remark, let us show that when D > 2, 
the exponent r is bigger than one. Let us assume that 
r < 1. Since th e sca ling function is integrable in zero, we 
can write Eq. ( 3.32 ) with a — 0, 



2 / f{x)dx 
lo 



f(x)f(y)K(x,y)dxdy (4.3) 



From the inequality, 

K(x,y) = (x% +y%){x~v +y~^) D ~ 1 > x x + y x (4.4) 
for D > 2, we see that Eq. ( |Q| ) leads to, 

f(x)dx> f(x)dx y X f(y)dy (4.5) 
o Jo Jo 



1G 



which, combined with Eq. ( |3.30 ), implies r > 1 in con- 
tradiction with our assumption. 



With the method briefly discussed in Sec. (BIB 3), it 



is easy to obtain exact bounds for r to control the results 
of the variational approximation we shall use. As a con- 
crete example, let us determine such bounds for D = 2 
and lo = 0.5. Since r < 1 + A (here, A = 0.75), Eq. 



(3.33) holds for a = A, for which we can numerically 
compute the minimum and maximum of A. From Eq. 
fl3.34|) , this leads to the inequality 1.5 < r < 1.607175. 
Thus, Eq. ( H| ) holds for 0.607175 < a < A, and we 
can compute new bounds for each a in this interval, and 
find the tightest bounds. The upper bound obtained for 
a = A cannot be improved since ^4(0) = 1 for a < A, 
hence 2 — (1 — a)/M a > 1 + a, but we obtain a better 
lower bound of 1.54 for a = 0.68. Table || presents such 
exact bounds for D = 2. 

For D = 1, the kernel is equal to x u + y u , correspond- 
ing to the kernel K\, with the notations of Eq. ( |3 . 3 1 ) , 
which was extensively studied in p6j . The exponent r is 
bigger than 1 for any lo > while t — when u) = 0. 
Since 1 < r < 1 + u> for lo > 0, we see that r — ► 1 when 
lo — > 0, hence r has a discontinuity at w = 0. 
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FIG. 8. The exponent r for the kernel (x" /D + y"/ ) - 1 
was computed with the variational approximation for lo > 
and 1 < D < 3. The theoretical w — > limit of r, to = 2—1/D 
is plotted on the Y-axis (squares). 



The variational approximation was used to study the 
behavior of r. The set of moments was chosen as dis- 
cussed in pq| . Results are shown on Fig. || for different 
values of D and lo > 0, while Fig. ^ shows the values 
of t for lo = 0. It appears that r has a discontinuity 
at to = not only for D = 1, but for £) > 1 as well: 
when lo — > + , r appears to have a limit To bigger than 
its value at lo = 0. Thus, the discontinuity which was 
rigorously shown to exist for D = 1, pertains for Z? > 1. 




It is difficult to accurately extract the value of To since 
the variational algorithm appears to be less accurate for 
small values of lo (for lo typically less than 0.1). However, 
To seems to be close to 2 — 1/D, which is the value of 1 + A 
at u) = 0. Actually, a heuristic argument, inspired from 
the discussion for the kernel in the large D (d > 1) 
limit, yields To = 2 — 1/D. 

Let fo(x) be the exact scaling function for lo = 0. From 
Eq. (EU0|), we get, 



TO 



lim 



+ OG 



(f(x)-f (x))x i+ — dx (4.6) 



and the limit in the right hand side of the latter equation 
must be strictly positive, although (/(x) — /o(^)) — * 
for any x > 0. How can this occur ? Since t > t w=0 
(for small lo), (f(x) — fo(x)) ~ c/x T when x — * 0, and 
c must vanish when cj — ► 0. Thus, the integral has an 
integrable singularity cx 1+ ^~ 1)/D ~ T . If r < 2 - 1/D 
(we know that To < 2—1/D from t < 1 + A) , the con- 
tribution of the singularity is wiped out by the vanishing 
of c, whereas, if To = 2 — 1/D, the integral is equivalent 
to c/(to +lo/D — t), and it has the finite limit To — r w= o 
provided that c vanishes as (to — t w= o)( t o + lo/D — t). 

Figure || plots the value of r and t = 2 — 1/D for 
a; = and 1 < D < 6. Both r and To have the limit 
2 when D — > oo which implies that the discontinuity in 
lo = vanishes at large D, as can be seen on the figure. 
The reason why r — > 2 is that when D — > oo, 



K(x,y) 



(xy)T+0(l/D) 



(4.7) 



and therefore the D — > oo limit of / = 2 D / is solution of 
Eq. (Al) with the kernel {xy)? which is a [i > kernel 
with exponent r = 2 (the same trick was used in p6| to 
study the d — > oo, d = AD limit for the i£"|> kernel). 



h 1.0 ' 




FIG. 9. Variational approximation for t when u> = 0, com- 
pared to its lo — > + limit to = 2 — 1/D. Both t and to tend 
to 2 when D — » oo. 
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For oj < 0, we have fi = ui/D < 0, and using the results 
of van Dongen and Ernst p5fl , we have, for a; — > 0, 

/(a:) - B^aT^ exp ^^^^ f +o ° x 1 - 1 ' 13 f(x)dx\ 

(4.8) 

where, _B, 7 and & = fim SS~ 2 ~ X Y are w dependent con- 
stants. These constants also depend on the definitions 
of Y(t) and S(t), but 7^2 when \i — ► 0. For a given 
definition, say 



5(t) =< s 2 > / < s >, Y(t) = S 2 /M ± , 



(4.9) 



van Dongen and Ernst showed that the scaling function 
crosses over to the \i = (polydispersed) case when 
uj — > 0, since the small a; asymptotics tends to B(0)a; _r , 
where r = 6 _1 (2 — J 00 x 1 ~ 1 / D f(x)dx) is precisely the 
w = polydispersity exponent (we had set b = 1 in Eq. 



(POD)- 

Consequently, we should observe this cross-over in nu- 
merical simulation. Moreover, for small, but finite u>, 
the critical x c below which f(x) is significantly depart- 
ing from the power law corresponds to /1 In x c of order 
one. Thus, it is reasonable to expect a scaling behavior 
when uj — ► 0~, 



f(x,w)=x c (u) T g(x/x c (uj)), 

with x c (oj) = exp(— c/w + o(l/u;)), <?(?/) 
and g(y) oc y~ r at large y. 



(4.10) 
at small ?/, 



w 


"7~m 


TM 


0.02 


1.020 


1.510 


0.2 


1.339 


1.588 


0.3 


1.472 


1.594 


0.4 


1.514 


1.601 


0.5 


1.540 


1.608 


0.6 


1.572 


1.614 


0.8 


1.623 


1.800 


0.9 


1.633 


1.900 



TABLE II. Exact upper and lower bounds 
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B. Numerical simulations 



Family and Meakin M and Meakin M found from sim- 
ulations that there was a polydispersity exponent for 
d = D = 2 and lo = 1/2, and interpreted the qualitative 
difference with the d < D regime as a gelling boundary 
effect. 

Indeed, for the actual heterogeneous growth model, it 
seems obvious that there is gelation for all d > D, while 
the mean-field theory finds that there is no gelation for 
A > 1, i.e. for d < D + 1 — lo. This tends to show 
that the mean-field theory fails at least for d > D. Nu- 
merical simulations performed for d = 2, to = — 1, and 
D = 1.7 > d + u — 1, show gelation at a finite time 
t g 0.12, which was seen to be nearly unaffected when 
doubling the mass of the sample, keeping the same ini- 
tial density and mass of the droplets, or reducing the 
time step by a factor 2, and thus seems to be well-defined 
in the continuous time and thermodynamic limit, corre- 
sponding to a genuine gelation transition. This is a con- 
firmation of the naive gelation criterion d > D. Actually, 
it is clear that the mean-field theory must break-down 
also for d = D. The reason is that the mean-field M\(t) 
diverges, while the actual M±(t) cannot diverge from a 
geometric constraint which, of course, is absent in the 
mean-field theory. Indeed, for d = D, Mi(t) is also pro- 
portional to the occupied area fraction, and is therefore 
bounded. This means that at large times, strong density- 
density correlations play an important role, and are not 
taken into account in the mean-field theory. 

In fact, the mean-field results may be qualitatively cor- 
rect in an intermediate regime between small times, and 
the asymptotic non mean-field scaling regime, and the 
scaling function may have a behavior qualitatively close 
from what expected for the mean-field. 

To check this, we performed simulations in d = 2, for 
various values of the growth exponent lo. In one step 
of simulation, all the droplets radii were increased of an 
amount 5r — r^St, then collisions were looked for and 
resolved. In most of the simulations, the time increase 
St was equal to 0.005. It was chosen small enough such 
that further reduction would not lead to significant mod- 
ification of the results. As can be intuitively understood, 
the number of droplets decreases much faster for d = D 
than for D > d. In the scaling regime, the number of 
droplets was cut by a factor of more than 1000, and we 
were obliged to start from a huge number of droplets 
(about 2.5 x 10 5 ) and to perform a large number of simu- 
lations to obtain acceptable statistics, without being able 
to reach very large times. Fig. [lO] shows a configuration 
obtained at t = 5.0 for to = — 1 from an initial configura- 
tion of 1024 2 droplets. It is striking that the distribution 
of droplets masses looks much broader than for D = 3 
(see Fig. ||). 




FIG. 10. A typical droplet configuration in the scaling 
regime of growth and coalescence with d = D — 2, obtained 
here for u — —1, from 262144 droplets of radius 0.75 in the 
initial condition on a 1024 x 1024 lattice. The picture repre- 
sents the whole system (with periodic boundary conditions) 
at t = 5.0 (S = 7693.9). The number of droplets has dropped 
to 287. 



The scaling form Eq. ( |3.17| ) was used with Y = tS 1+p 
to obtain convincing data collapse, as shown on Fig. |i~l| . 




0.05 



s/S(t) 



FIG. 11. Scaled mass distributions at three stages in the 
simulation of heterogeneous growth with u) = — 1. These re- 
sults were obtained from 96 simulations. In each simulation, 
262144 droplets of radius 0.75 were initially randomly placed 
on a 1024 x 1024 lattice (without overlap). 



Fig. [12] plots the scaling functions for various values 
of lo. The results are consistent with a transition from 
a small x diverging scaling function, for lo > 0, to a 
small x vanishing scaling function for negative value of 
lo. For the considered values of negative lo, the scaling 
function, as visible on Fig. although vanishing when 
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x — > is quite broad, with a maximum at a value x sig- 
nificantly smaller than 1. When u> — > Q — , we observe a 
cross-over to the u> = power law, and the position of 
the maximum of the scaling function rapidly tends to zero 
when lo — ► 0, consistently with the discussion around Eq. 
(4.10). Moreover, the exponent extracted from the nu- 



merics is about 1.2, which compares well with r = 1.108 
from mean-field. However, the r exponents for u> = 
and u> = 0.5, do not seem to be significantly different, in 
contrast with the quite large discontinuity in mean-field. 



4.0 



3.0 - 



2.0 - 



1.0 



slope 0.5 



As far as S(t) is concerned, it is difficult to be positive 
due to strong numerical limitations. Figure [l3| shows 
the evolution of S(t) for uj = —3. The excellent data 
collapse shown in Fig. |lj corresponds to t = 65.0 and 
t = 77.5. At these times, S(t) is seen to grow much faster 
than t £> /( 1 ~") = -y/t, which is consistent with the mean- 
field logarithmic enhancement. This may explain why 
the mean-field description is qualitatively correct in the 
scaling regime observed in simulations. A strongly non 
mean-field scaling with properties closer to the d < D 
case may be observed at times unreachable to our simu- 
lations. 



o.o 



-1.0 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

iog(t) 

FIG. 13. For u = -3, S(t) grows much faster than t D/ ^~ u) 
(slope 0.5) even in the scaling regime (t > 50) corresponding, 
a behavior which may be related to logarithmic enhancement 
of S in mean-field. 



V. CONCLUSION 



o 
o 







o h 



_CD=0 



(0=0.5 




m=-3 



-5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 

log(s/S(t)) 

FIG. 12. Small x — s/S behavior of the scaled mass distri- 
butions obtained in numerical simulations for different values 
of the growth exponent u. 



In this article, we have extensively studied a gener- 
alized Smoluchowski equation corresponding to aggre- 
gation processes for which particles (or clusters) grow 
between collisions, with s = As@, and small particles 
(monomers) are injected. 

A physical motivation for this work is droplets nucle- 
ation and we have directly derived generalized Smolu- 
chowski equation and found the collision kernel for two 
models, respectively describing homogeneous and hetero- 
geneous nucleation. 

For a generic kernel, with parameters A, /i, we have 
shown that the gelation criterion was max(A, (3) > 1. We 
have devoted much time to the study of the equation 
without injection, for which we have provided two ex- 
act solutions. The scaling properties for a generic kernel 
are seen to be strongly affected by the exogenous growth 
term, and depend on j3, A, and (i. For A > 0, however, the 
scaling is the same as for the standard equation. For the 
interesting case A = 0, the behavior of the typical mass 
S(t) is modified, but the scaling function is unchanged. 
For A < P, the scaling function is qualitatively different, 
and vanishes at a finite xq > 0. 

We have also studied the case of a constant injection 
rate of monomers. The distribution reaches an asymp- 
totic steady state with a power law tail N QO (s) cx s~ r , 
and we find that t depends on (3 and A. We have shown 
that this steady state can be seen as a sol-gel critical 
state, breaking the mass conservation by collisions, which 
we think clarifies our physical understanding of these 
power law asymptotic states. 
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We have payed special attention to the case of a con- 
stant mass injection rate and A = 2/3—1, related to 
homogeneous nucleation. This corresponds to a time- 
dependent, self-consistent injection rate I(t) oc c—Mp(t). 
We have shown that I(t) vanishes at large time, in agree- 
ment with the droplets deposition and coalescence model. 
For droplets deposition, Mg is proportional to the surface 
coverage, and the vanishing of the injection rate corre- 
sponds to the saturation of the coverage to 1. Thus, our 
self-consistent Smoluchowski's equation recovers a merely 
geometrical constraint, which is quite nontrivial. 

As far as scaling is concerned, we have found consis- 
tent results for fi < (3 — 1, with nontrivial polydisper- 
sity exponents for fi = (3—1, recovering 9 = 1 + [3 
and z — 1/(1-/3) as for the droplets deposition and 
coalescence model. However, for fi > [3 — 1 kernels, we 
could not find a consistent scaling, and there might be 
no scaling solution with a constant mass injection rate. 
The mean-field kernel for droplets deposition and coa- 
lescence has ji > /3 — 1, and we argued that taking into 
account excluded volume pair correlations may be essen- 
tial to obtain a consistent description by a kinetic equa- 
tion, including a nontrivial polydispersity exponent if /x 
is switched to j3 — 1. It is quite difficult however to study 
these correlations either numerically or analytically. 

Finally, we have applied these results to droplets 
growth and coalescence with d = D. We have shown that 
Smoluchowski's approach accounts for the qualitative dif- 
ference in the scaling function with the d < D case. We 
have computed nontrivial polydispersity exponents oc- 
curring for uj > 0, and described the cross-over from 
monodispersity to polydispersity occurring for u> — ► CP. 
We have compared these theoretical results with numeri- 
cal simulations, with good agreement, despite mean-field 
limitations which we discussed. 

As a conclusion, we would like to point out that one 
of the main reasons why people have become increas- 
ingly interested in breath figures is that it is an exam- 
ple of a geometrically constrained growth process, where 
diffusion plays a minor role, in contrast with diffusion 
limited cluster-cluster aggregation p8|-f4(| or Brownian 
coalescence of droplets Q. Therefore, one could doubt 
that neglecting density-density correlations may have no 
dramatic consequences. Indeed, for homogeneous nucle- 
ation, we have seen that pair correlations may be cru- 
cial to find a correct scaling function, and we also found 
an infinite upper critical dimension for heterogeneous 
d = D nucleation. However, we have shown that Smolu- 
chowski's equation in an extended form could be suc- 
cessfully used to describe heterogeneous growth, and also 
gives very interesting insights into homogeneous growth, 
which was not a priori obvious. 

We are very grateful to P.L. Krapivsky for helpful cor- 
respondence. 



APPENDIX A: SCALING THEORY WITHOUT 
INJECTION 

In this appendix, we shall give a detailed demonstra- 
tion of the scaling results given in Sec. Ill B 2 for gen- 
eralized Smoluchowski's equation without injection of 
monomers. We make the assumption that starting from 
a monodispersed distribution of droplets, the late time 



solution of Eq. (3.1) has the scaling form of Eq. (3.17). 
Note that, although it is quite clear from its homogene- 
ity behavior that Smoluchowski's equation admits scaling 
solutions, it has never been mathematically proved that 
these solutions are approached at large times, except in 
the few cases for which we can obtain the exact solution. 
However, scaling is commonly observed experimentally 
and numerically for aggregation models, as well as in nu- 
merical solutions of Smoluchowski's equation |3(| (when 
possible), making this assumption very reasonable. 

Another important remark is that Smoluchowski's 
equation may generally speaking admit different classes 
of scaling solutions, with different asymptotic behav- 
ior, and that one peculiar class will be selected from 
the initial condition. This must be remembered in the 
discussion. In our demonstration, we shall always use 
the fact that S(t) cannot be negligible compared to 
So(t) ~ t 1 ^ 1 - ^), the lower cut-off, from the discussion in 



Sec. Ill B 2 . We shall also implicitly assume that Y(t), 
S(t) and all the moments of N(s, t) are asymptotically 
regular convex or concave functions. 



1. Within the nongelling domain 

To start with, let us study the case A < 1 and j3 < 1, 
i.e. nongelling systems which are not on the gelling 
boundary. Our discussion is based on the fact that ei- 
ther S(t) > s (t) or S(t) oc s (t). We shall study the 
implications of both possibilities. 



a. Case S 3> so 

Let us assume that S(t) ^> s (t), or equivalcn tly t hat 
S ,/3_1 < S/S. We see t hat the growth term fl3.25| ) is 
much smaller than ( |3 . 24 ) in the scaling limit, and the 
scaling equation is, 



bxf'(x) + af(x) = f{x) J +o ° f{x\)K(x, X\)dx\ 

~ klo f( x i)f( x - xi)K(xi,x - xi)dxi (Al) 



-l-A 



where b = lim SS~ Y and a = lim YS~'~" are posi- 
tive, and possibly zero or infinite. 

For finite a and b, this equation is very close from stan- 



dard Smoluchowski equation Eq. (1.1) with the same 
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kernel K, which corresponds to a — 2b and < a < +00 
and was well studied in the literature p^p6|. The 
polydispersity exponent r, if any, has the upper bound 



t < 1 + A < 2 (read Sec. IIIB3). As a consequence, from 
Eq. (|3.2lD we have, 



M x {t) oc S 2 /Y 



(A2) 



From Eq. fl3.27| ), we see that Mi(t) is nondecreasing. 
Thus, either Mi(t) tends to a finite limit, or it goes to 
infinity. 

(i) If Mi(t) tends to a finite limit, then, necessarily 
from Eq. (|A|), Y cx S 2 . The scaling of ( |3~2^ ) with 
( |3.26 ) requires S cx S x , hence, 



S(t) oc t~ 



(A3) 



To be consistent with our assumption that S(t) 3> 
fVC 1 -/ 3 ^ we must have A > 0. Besides, from Eq. ( [3.27[ ), 
a necessary condition for Mi(t) to have a large t finite 
limit is that Mp(t) must be an integrable function. If 
r < 1 + 0, M scales as S 1+l3 /Y, i.e. M cx t^- 1 )/^-*) 
and is integrable since A > /3. 

If t > 1 + 0, we have, from Eq. (|3.22| ), 



M^ a 5(t) T - 2 s (t) 



1+/3-T 



< s (i) /3 " 1 cx t 



(A4) 



since r < 2 and S'(i) » So(t). Therefore, Mg, being 
equivalent to a power law <C is integrable. 
If r = 1 + /?, from Eq. ( ^23| ), 



(A5) 



Mg cx SXO^ 1 ln(5/s ) oc lni 

which is integrable when A > 0. 

(ii) Now, let us consider the case when M\(t) diverges 
at large time. From Eqs. (|3~27| ), (|A§), ( fj| ) and ( ^23| ), 
we obtain, 



MiS"- 1 , ifr<l + /3 

Mi cx < A/ 1 5 T - 2 s (i) 1+ ' 3 - r , if r > 1 + /? (A6) 
Mi Saints'/so), ifr = l + /3 



Anyway, in the three cases, S 3> so implies that Afi <C 

Mi 



MiSg 1 , hence, 



Mi < 



(A7) 



Therefore Mi(t) < t Q for any a > 0. From Eq. (|A§) 
and the fact that SyS* is at least of order 1/t since 
S(t) > tV(i-« ; Eq. @ requires that, 



S Y 



(A8) 



Thus, the scaling condition between ( p. 24 ) and ( 3.26| ) is 
simply, 



cl+A 



(A9) 



which implies that S 1 A is dominated by a power law. 
Combined with the fact that S(t) > this re- 

quires that 



and 



S 1 



Mi oc t^S 1 '* > t- 



(A10) 



(All) 



thus, from Eq. (A7), we m ust have A > 0. Now, com- 
bining Eqs. (All) and (|A6|), we see that, 



Mi oc Mf^i-^ln^ wrMi))" 3 



(A12) 



(oc3 = lifr = l+/3, otherwise 03 = 0). 

Since Mi — > 00, the right hand side of the latter equa- 
tion must be non integrable, and as Mi is much smaller 
than any positive power of t, this implies that a 2 < 1 . 



a 2 




if A < 
if t > 1 + (3 



(A13) 



Since r< 1 + A, a 2 > 1 if A > /3. Therefore we must 
have X < 0. However we already found that X> 0, thus 
A must be equal to /3. 

As a consequence r is never bigger than 1 + = 1 + A. 
and we can distinguish between fi > kernels, for which 
t = 1 + X = 1 + and fi < kernels for which there is 
no polydispersity exponent or r < 1 + 0. 

Let us star t wi th fi < kernels. As r < 1 + /?, and 
A = 0, Eq. ( |A12|) is reduced to Mi oc 1/t, which leads 
to, 



Mi(t) oc Int. 

For a, fi > kernel, r = 1 + A = 1 
leads to, 

Mi oc (lnMi)/£, 

and it is easily seen that, 

Mi(<) oc (mi)ln(lnt). 

In both cases, we have, 

S(t) oc (tMx)T^, 
F(i) oc S 2 /Mi. 



(A14) 

/?, and Eq. (|Al|) 
(A15) 

(A16) 



(A17) 
(A18) 



Thus, the initial assumption that S(t) ^> so(t) implies 
that A > and that the scaling equation is Eq. ( Al) 
with a = 2b (since 2S/S ~ Y/Y) and < a < +00, 
and is the same as for standard Smoluchowski's equation 
with the same kernel. 
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b. Case S oc so 



Conversely, let us assume that S(t) oc t 1 ^ 1 ^^. If 
Y/Y ;§> S/S oc 1/t, Y increases faster than any power 
law, the growt h ter m is still negligible at large time, 
and Eqs. ( 3.24 ) and ( 3.26 ) are of the same order, hence 
S 1+x /Y oc Y/Y, which is contradictory, for S 1+x /Y van- 
ishes faster than any power law while Y/Y ^> 1/t. 

Thus, Y/Y = 0(S/S) and both terms in the right 
hand side of Smoluchowski's equation are of the same or- 
der at large time, as S/S oc S*^" 1 . Besides, both terms 
must scale as the collision term, otherwise the obtained 
scaling equation has no physical solution vanishing below 
a finite argument xq > 0. Thus, ( 3.25| ) and ( 3.26| ) must 
be of the same order, which yields: 



Y 



-1 5 1 + A 

OC 



Y 2 



(A19) 



and Y(t) oc S'(£) 2+A ~' 3 . The fact that the scaling function 
vanishes at a finite x^ > ensures that, M\{t) oc S 2 /Y, 
hence, 



Mi(<) oc S{t) 



0-X 



(A20) 



Since M\{t) is non decreasing, we must have A < /3. How- 
ever, if A = (3, Eq. (13.3!) yields 



(A21) 



Mi oc 1 /t, 



hence M\(t) oc hit, which is in contradiction with Eq. 
flA2C| ). Therefore one must have A < (3. 
The scaling equation has the form, 



b [0f(x) + xf'(x)} - a(x^f(x)Y = f(x) / f(xi)K(x, Xl )dxi 



+ 00 



f(xi)f(x - x 1 )K(x 1 ,x - xij dx 1 



(A22) 



c. Conclusion 



2. A = 1 and (3 < 1 



Since S(t) > so(t), the collection of the two cases we 
examined above leads to the conclusion that, if A < 1 
and p < 1, there are three main regimes of scaling, in 
agreement with the qualitative discussion in Sec IIIB2. 

If (3 > A, S(t) scales as S g (t) oc iVU-fl, Y(t) oc S(t)», 
with 9 = 2 + A — (3, and there is no polydispersity ex- 
ponent since the scaling function is zero below a finite 

Xq. 

If [3 < A, S(t) scales as S c (t) oc t 1 ^ 1 ^^, the mass is 
asymptotically conserved, i.e. Mi(t) tends to a constant, 
and 9 = 2. There can be a polydispersity exponent, 
which is t he same as for standard Smoluchowski's equa- 
tion Eq. ([□]) with the same kernel. 

Eventually, in the marginal case when \ = [3, the scal- 
ing of S(t) depends on the kernel not only through its 
homogeneity A, but also through its /i exponent. As in 
the (3 < A case, the scaling equation is the same as for 
Smoluchowski's equation with the same kernel. For /i < 
kernels, the mass in the system Mi(t) oc lni, while for 
H > kernels, with t = 1 + A, Mi(t) oc t(lni) ln(lnt). In 
both cases, S oc (iM 1 ) 1 /( 1 ~' 3 ) and Y oc S 2 /M x . 

These scaling results can be compared to the case K = 
1,(3 = which we solved exactly: we found that the con- 
ventional scaling breaks down and that, with the proper 
scaling form, the scaling function is the same as for the 
exactly solvable standard Smoluchowski's equation with- 
out the growth term, that S(t) oc (tint), M\(t) oc lnt 
and Y oc t 2 lni, just as predicted by the scaling theory. 



For A = 1 and [3 < 1, it is possible to follow the same 
line of reasoning, with a few modifications. In this case, 
one has to distinguish between fi > and /x > (this 
is also true for standard Smoluchowski's equation with 
A = 1 |25|]), since for /i > 0, we find r = 2 and the scal- 
ing of Mi has an extra ln(S'/so)- It is found that Mi 
is asymptotically conserved, that S(t) 3> so(t), and that 
the scaling equation is Eq. (Al) with a = 2b < +oo. For 
/i > 0, one has S oc S/(lnS), which leads to, 



whereas if fj, < 0, 



S(t) oc e bVi , 



S(t) oc 



(A23) 



(A24) 



where b cannot be derived from the scaling theory. 



3. (3 = 1 and A < 1 

In this case So(t) oc e', and the discussion is quite dif- 
ferent. From Eq. (3.27), we see that: 



Mi(i) = Mi(0)e* 



(A25) 



Since S(t) > so(t) oc e*, we have in the large time limit 
S/S > 1. Let us assume that S/S ^> 1, i.e. that S(t) is 
bigger than any exponential function e at , which entails 
that (3.24) is much bigger than (3.25) (which scales as 
l/Y since (3 = 1). From Eqs. (|3.21|), ([3~22|) and ( p3| ), 
it is clear that, 
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S 2 /Y = 0{M 1 {t)) = 0{e t ). 



(A26) 



4. A = 1 and fi = 1 



Consequently, if ( |3 24j ) scales as (|3.26|) , then S/S x = 
0(e*) and y/yU+Wa = o{e t / 2 ). Since A < 1, these two 
relations are in contradiction with the assumption that 
5 is much bigger than any exponential function (whi ch 
implies the same property for Y , through Eq. ( A26 )). 
We see that if 5/5 > 1 , (|3.24| ) is the leading term in 
the scaling limit, and the scaling function is a pure power 
law f(x) = cx- T with r = Hm(Y /Y)(S/S). One must 
have r > 2 such that the total mass in the system be 
finite at finite time in the scaling regime. Making use of 
Eq. fl3.22| ), we find that n{t) oc Mi(t)/s (t) would tend 
to a finite value > 0, which is unphysical. 

Therefore 5/5 is of order 1. From arguments very 
similar to those we just used, it is easily seen that Y/Y 
cannot be much bigger than 1. Thus, (3.24) scales as 
( p5| ). If Y/Y -> 1 and S/S -> 1, the left hand side 
of Eq. (3.1) vanishes and we have to take into account 
the subleading terms in the scaling limit. This occurs for 
A = as will be seen below. 

If t he left hand side does not vanish, the scaling with 
( 3.26 ) lead s to 5 1+A oc Y, and the scaling equation is once 
again Eq. @ , but now with b/a = (1-S/S)/(1-Y/Y). 
Consequently, the polydispcrsity exponent, if any, is less 
than 1 + A, and Eq. (|a|) holds, leading to S 2 cx Ye 1 . 
Since S 1+x oc Y, we have, 



5(t) oc e~ 
Y(t) oc S 2 e- t 



(A27) 
(A28) 



which excludes A < since S(t) > so(t) — e*, and also 
A = for which S/S -> 1 and Y/Y -> 1. Note that in 
the A > case, we find a = 2b, and once again the scal- 
ing equation is the same as for standard Smoluchowski's 
equation. 

Indeed, for the exactly solvable case K = 1, ft = 1, 
which corresponds to A = 0, we found that S(t) oc te* 
and Y(t) oc S 2 /e\ thus S/S -> 1 and Y/Y -> 1. Thus, 
to treat the A = case, we shall write S(t) = X(t)Mi(t), 
with X/X < 1, and we have Y oc S 2 /Ah = XS. The 
right hand side of Eq. (3.1) scales as, 



- 1(2/(3;) 
while the left hand side scales as 
£(•■)«£(..) 



(A29) 



(A30) 



which leads to X(t) oc t, recovering the exact result for 
K = 1. Once again the scaling function is Eq. (Al) with 
a — 2b < +oo. The polydispersity exponent r is strictly 
less than 2, which justifies a posteriori that Y oc S 2 /Mi 
(it is possible to show that assuming r > 2 leads to a 
contradiction). 

However, for A < 0, we were unable to find a consistent 
scaling. 



In this case, we still have M±(t) oc e* and so(t) ~ e*, 
but it is easily seen with the same kind of arguments as 
above, that o ne mu st have S/S ^> 1. Thus the exog enou s 
grow th te rm ( 3.25| ) is negligible, and the scaling of ( 3.24 ) 
and flH|) yields, S/S oc S 2 /Y. 

For fi < kernels, one has Mi oc S 2 /Y and we obtain 
S/S oc e*, leading to, 



S(t) 



(A31) 



For n > 0, we have r = 2 and M x oc S 2 ln(5/e*)/Y, 
leading to 5/5 oc e*/ln(5), and, 



5(t) oc e 



(A32) 



In these expressions b is an unknown positive constant. 
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